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Abstract 

In this paper we extend the estimate of the value of the lattice spacing a in 
units of the ro scale at values of the bare coupling larger than those avail- 
able. By using results from the computation of the renormalised coupling in 
the Schrodinger functional formalism we find that from (5 ~ 7 onward the 
behaviour predicted by asymptotic freedom at tree loop describes very well 
the data if a value of roA slightly lower then the latest available is used. We 
also show that, by sticking to the current value, an effective four loop term 
can describe as well the data. The systematic relative error on the lattice 
spacing induced by the choice of the procedure is between 1% and 3% for 
6.92 < < 8.5. 



1 Introduction 



In order to obtain physical predictions from a lattice QCD simulation, one 
has to spend as many "experimental" input as are the free parameters of 
the theory. This situation is not peculiar of lattice regularisation, but is a 
common fact to every regularisation/renormalisation procedure. In the pure 
Yang-Mills case the only free parameter of the theory is the bare coupling 
go, and the only dimensionful quantity in a simulation, in the idealized case 
of a lattice of infinite extent, is the lattice spacing a. Asymptotic freedom in 
the high energy regime predicts the functional dependence of a from go once 
Al is known. If the value of the bare coupling gets too large to rely on a 
perturbative expansion, one has to resort to a non-perturbative definition of 
the physical scale of the simulation. 

One way to perform this task, in the case of a lattice extension large 
enough to safely account the dynamics of the light quarks, is to use the ex- 
perimental inputs coming from the spectroscopy of the strange mesons, as 
explained in refs. pi pi H. Another way to set the scale, widely used since its 
proposal in is based on the calculation of the force between static colour 
sources; the physical input in this case is a length scale, ro — 0.5 fm. 

In H and then in || the quantity a/ro has been computed, in the quenched 
QCD framework, for a range of j3 = Q/g^ going from (3 m \ n = 5.7 to /3 max = 6.92 
and a parametric description of ln(a/ro) as a function of g\ is provided that 
describes the data with an accuracy better than 1% in the whole range. In the 
quenched approximation the absolute value of the scale is affected by a sys- 
tematic error, induced by the quenching, of the order of 10%; nevertheless the 
precise result of ref. ||] has been widely used in literature, at least to set the 
ratio between different scales and to perform continuum extrapolations. The 
use of the parametrization of ref. || outside its validity range (5.7 < (3 < 6.92) 
is highly questionable; on the other hand, one can easily encounter situations 
(see for example ref. ||) where a precise knowledge of a/ro f° r values of f3 
greater than 6.92 is required. 
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In this paper we will show how to obtain an accurate estimate of a/ro 
for (5 > 6.92, using existing numerical simulations and the behaviour expected 
from asymptotic freedom: we use the non-perturbative results coming from 
the computation of the Schrodinger functional renormalised coupling flO| , 

1 l| and the parametrization dictated by Renormalisation Group improved 
perturbation theory. 

2 Non— perturbative approach: the running coupling 

In a series of publications (see for example refs. []|, |l0|, [ll], |l2|) it has been 
shown how to define and compute a renormalized coupling, whose running 
with the energy scale can be followed, in the continuum limit, over a very 
large range, making able to connect the high energy perturbative region with 
the low energy non-perturbative one. 

The calculations has been done using a recursive finite-size technique 
in the Schrodinger functional framework where the energy scale is given by 
fjL = 1/L. The size of L m3iX , defined implicitly by the relation 5g F (L max ) = 3.48, 
has been computed in ref. || in terms of ro, with the result 

WAo = 0.738(16) ; (2.1) 



in the same paper, combining this result with those of ref. [11], the authors 
quote 



A^ = 0.586(48)/r , (2.2) 



MS 

where the superscript (0) remind us that we are dealing with the quenched 



approximation. Equation (|2.1|), together with the tab. (6) of ref. [11|, will 
be the handles to extend our non-perturbative knowledge of a/ro in a wider 
range of (3 values. 

We start by considering a length scale L u , implicitly defined by the rela- 
tion g 2 (L M ) = u: by construction, this scale is smaller than L max by a known 
factor s u . We can use this fact in order to express the scale L u in terms of rg: 
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^ = x ; (2.3) 
if we consider a particular discretisation of L u , i.e. L u /a = N, we can write 

a/r = x/N . (2.4) 



From table (6) of ref. [11| we can read off a whole bunch of pairs (L u /a, (3) 
and, using eq. (|2.4j) , we can compute a/ro corresponding to all these values of 

Given the finite extent used in these calculations we must, of course, 
assume that lattice artifacts do not play a prominent role. Whether this is a 
good assumption can be checked by using the /3-value corresponding to the 



coupling u = 2.77 in tab. (6) of ref. [11|, that lies within the range in which 
a/ro has been directly computed. In this case the scale factor s u is almost 
exactly equal to 3/2. 

The results can be seen in table ([J): three of the values of a/ro computed 
through eq. ( |2.4j ), including the points computed at u = 3.48, are at j3 < 
6.92, and nicely fall, within the error, which comes entirely from the error on 
Lraax/i"o, over the parametrization of ref. |]], even if the points correspond to 
lattices with a relative poor discretisation, N = 8 and N = 12. 



Indeed, the results of ref. [11] for the continuum limit of the step scaling 
function also indicate that lattice artifacts are very small already at N = 8 
and N = 12, strongly supporting our assumption. 



p 


eq.(2.6) ref. [6] 


eq. (2.4) 


it 


Lu/a 


6.4527 


-2.345 


-2.383(22) 


3.48 


8 


6.7750 


-2.758 


-2.790(22) 


3.48 


12 


6.7860 


-2.773 


-2.789(22) 


2.77 


8 



Table 1: Check of eq. (2.4) against non-perturbative result of ref. [6]. The accuracy on 
the parametrization of ref. [6], in this f3 range, is about 1%. The error on the results of eq. 
(2.4) comes entirely from the error on I/ ma x/ro. 
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3 Perturbative approach 



In ref. [13] it is argued that the three-loop perturbative expansion in the qq— 
scheme may be trusted up to a q q ~ 0.3. We test how perturbation theory 
can help us in safely extending the region where a/ro has to be considered 
as known with good precision, contrary to the general belief that lattice bare 
perturbation theory is well-behaving only in the extremely huge cutoff region. 

We recall some well known results. In a generic renormalisation scheme, 
S, the perturbative expansion of the /3-function reads: 

^9s = P(9s) = -tslbo + hf s + b[ S) t S + •••} (3.1) 

where fro an d b\ are the well-known universal terms, 

12 1 38 

(S) 

while, from b 2 on, the coefficients of the /3-function start to be scheme- 



dependent. The b 2 coefficient is known in the MS scheme |17| : 



AMS) 1 ,77139 5033 325 7 . , . 

Since the relation between and the bare lattice coupling go is known up to 
two loops Jl4], |l5[| , some straightforward algebra leads us to the b\ coefficient 
of the lattice /3-function, that we quote in the Nj = case only: 

b\ = bf* + Mi(1)/(4tt) + {d\(l) - d 2 (l))6 /(4vr) 2 = -0.0015998323314(3) 

(3.4) 

where d\ and c?2 can be found in ref. |l4| . 

In order to integrate the renormalisation group equations, we need the 
integration constant, usually defined trough the renormalisation QCD scale 
Ag, which can be related to the running coupling g s (/i) through the exact 
relation 

A S = „teir"«>e-/«l> exp {_ £ d4 _L + ^ _ |L] }. ,3,, 
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where the subscript S indicates a generic scheme. The previous expression 
obviously becomes a perturbative one upon the insertion of the perturbative 
formula ( |3.1| ) for the /3-function. In the case of lattice bare perturbation 
theory, we can recast the previous equation in the following form: 

In(o/r ) = - ln(A L r ) - i ln OoSo) - tJ~^ ~ Hdo) , (3-6) 
where I (go) is the integral appearing in the exponential of eq. (|3.5|), By using 



eq. ( |2.2D and the result from ref. [16] 



.,0, 

we get 



A^/Al = 28.80934(1) (3.7) 



A L r = 0.0203(17) (3.8) 

and, with this input, eq. (3.6) becomes a perturbative tool to set the scale a 
in terms of ro- 

As can be seen in fig. (|l|), the three-loops formula describes quite well the 
behaviour of non-perturbative data at high [5 and nicely superimpose to the 
points computed through eq. (2.4) if the value of A is lowered by about 1 a 
with respect to the value quoted in ref. || (the upper bound of the three-loops 
band in the figure corresponds to A^L = 0.538/ro). 

If we do not stick to the value A^L = 0.586/ro but let the three- loops 
fit decide for its value, we obtain the results listed in table [2|. The value 
of A^ obtained through this procedure is remarkably constant against the 
variation of the number of points included in the fit. Whether this fact is really 
suggestive of a lower value of the A parameter in quenched approximation 
could be only clarified with further non-perturbative investigations at smaller 
values of the bare coupling. 

Fig. (Q) indicates that, at [3 ~ 7, three-loops perturbation theory is not 
far from an accurate description of non-perturbative data even when we keep 
A^ = 0.586/ro. We can thus try to determine an effective four-loops term 
by means of a fit. A two parameters fit (A and 63) gives a value for A^ro 
totally compatible, in the error, with 0.586, and this make us confident that 
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Figure 1: ln(a/ro) as a function of gl. 



reducing the number of parameters by sticking to the known value of Al^ro 
will give a reliable estimate of 63. 

In table (|3|) we show the results obtained with a one-parameter fit, namely 
6| ff itself, by keeping fixed A^ = 0.586/Vo. We note that absolute value of bf 
agrees with an approximate geometric growth of the coefficients (47r) 2 ^ n+1 ^6 n . 



n A^r ln(a/r )(/3 = 7.5) ln(q/r )(/3 = 8.5) 

3 0.547 -3.640 -4.801 

4 0.546 -3.638 -4.800 

5 0.540 -3.627 -4.789 

6 0.535 -3.618 -4.779 



Table 2: Results for the 3-loop fit to non-perturbative data at large (3, including the last 
four points obtained through eq. (2.4); n is the number of points included in the fit (most 
perturbative ones). 



n 



bf ln(a/r )(/3 = 7.5) ln(a/r )(/3 = 8.5) 



3 
4 
5 
6 



-0.0022 
-0.0022 
-0.0023 
-0.0025 



-3.643 
-3.643 
-3.638 
-3.631 



-4.819 
-4.818 
-4.815 
-4.810 



Table 3: Results for the 4-loop fit to non-perturbative data at large (3, including the last 
four points obtained through eq. (2.4), keeping fixed A^lro = 0.586; n is the number of 
points included in the fit (most perturbative ones). 

Summarizing, it is safe to assume that beyond (3 = 6.92 the quantity 
ln(a/Vo) is well described by an effective four-loops perturbative curve, if not 
already by the three-loops expression. We give as a final result the value 
obtained with n = 6 and A^r fixed to the value 0.586 {bf = -0.0025(3): 
the x 2 /n.d.f. of the fit is 0.9), that we chose in order to obtain the lattice scale 
a in terms of ro for j3 > 6.92: the fit is shown in fig. (|l|). The accuracy of 
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this description can be estimated by considering the variation of this quantity 
when changes by 1 a. Of course the value of will be readjusted by 
the fit. We take as a reference points the values of ln(a/ro) at (3 = 7.5 and 
(3 = 8.5, where we obtain respectively (assuming tq = 0.5 fm) 

P = 7.5, ln(a/r ) = -3.63(2) a = 0.0133(3) fm, 

(3 = 8.5, ln(a/r ) = -4.81(3) a = 0.0041(1) fm (3.9) 

We can therefore state that, no matter which solution we adopt, lower A and 
3-loops or standard A and 4-loops, the uncertainty on a/rQ for 6.92 < j3 < 7.5 
can be bounded by 2%, raising to 3% at f3 = 8.5 and asymptotically reaching 
the value of the uncertainty on A of 8%. 
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